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Computer simulation of coherent interaction of charged particles and photons with 

crystalline solids at high energies 

Armen ApyarQ 

Northwestern University, Department of Physics and Astronomy, 2145 Sheridan Road, Evanston, IL 60208, USA 

Monte Carlo simulation code has been developed and tested for studying the passage of charged 
particle beams and radiation through the crystalline matter at energies from tens of MeV up to 
hundreds of GeV. The developed Monte Carlo code simulates electron, positron and photon shower 
in single crystals and amorphous media. The Monte Carlo code tracks all the generations of charged 
particles and photons through the aligned crystal by taking into account the parameters of incoming 
beam, multiple scattering, energy loss, emission angles, transverse dimension of beams, and linear 
polarization of produced photons. 

The simulation results are compared with the CERN-NA-59 experimental data. The realistic 
descriptions of the electron and photon beams and the physical processes within the silicon and 
germanium single crystals have been implemented. 



PACS numbers: 41.60.-m, 25.20.Dc, 24.10.Lx, 87.18.Bb 
Keywords: Monte Carlo simulations, Coherent Bremsstrahlung 



Single Crystal, Electromagnetic shower 



I. INTRODUCTION 

In the last decade, the electromagnetic interaction of 
charged particles and photons with crystalline and amor- 
phous media are intensively investigated both theoreti- 
cally and experimentally. Special attention was given to 
the investigation of interaction of charged particles and 
photons with crystalline media. An electron or photon 
impinging on a crystal will interact coherently with the 
atoms in aligned crystal axes or planes. If the Laue con- 
dition is satisfied, the coherent bremsstrahlung (CB) or 
coherent pair production (CPP) phenomena are mani- 
fested The essential characteristics of the phenom- 
ena are quasi-monochromatic spectrum, high intensity 
and linear polarization degree of radiation in coherent 
maximum. The intensity of the coherent bremsstrahlung 
in aligned crystals is few ten times greater than in amor- 
phous media at the certain energy region. The same 
characteristics are seen in CPP by photons in the single 
crystals. The processes of CB and CPP in single crys- 
tals are well investigated and understood both theoreti- 
cally Q and experimentally 0, Q . 

In this work we describe briefly the physical processes 
involved, the simulation model and the results of the 
Monte Carlo simulations for high energy particles travers- 
ing an aligned single crystals. 

It is well known that the dominant energy loss mecha- 
nism for high energy electrons and positrons is the pro- 
duction of electromagnetic radiation, i.e. bremsstrahlung 
for motion through the matter. The high energy photons 
become absorbed mainly due to the e + e~ pair produc- 
tion in a matter. This is true for crystals as well. In 
our simulations we ignored the processes such as Comp- 
ton scattering, photoeffect, energy loss in ionization, nu- 
clear processes due to smallness of their cross section and 



negligible contribution in total cross section. Electron 
or photon beams penetrating the crystal will create an 
electromagnetic shower (EMS) via the CB, CPP or inco- 
herent bremsstrahlung (ICB) and pair production (IPP) 
depending on the orientation of crystal axes and planes 
relative to incident particle momentum. 

There are broad experimental and theoretical inves- 
tigations devoted to the EMS development in various 
amorphous media and energies. Many general - purpose 
Monte Carlo simulation packages (see Ref. [i| and refer- 
ences cited therein) for transport of particles and radi- 
ation through the amorphous media are developed and 
have been successfully used in the last decades. 

The electromagnetic interaction processes in aligned 
crystals are more complicated than in amorphous media. 
The cross section of interaction strongly depends on the 
crystal type, its orientations as well as on the energy, 
angular distribution and polarization of initial photon or 
electron beams. 

The EMS in single crystals was mainly considered in 
channeling regime, when electrons and photons penetrate 
the crystal along a direction close to the one of main 
crystallographic axes @] • We consider EMS in single 
crystals oriented in CB mode. All the above mentioned 
peculiarities of coherent processes are carefully taken into 
account in Monte Carlo computer code. 



II. THEORETICAL BACKGROUND 

Differential cross section of coherent radiation in a 
crystal is composed of two terms [l[ coherent and in- 
coherent bremsstrahlung 



da = d(jj 



da r 



(1) 
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The first term corresponds to the incoherent cross sec- 
tion (including radiation in the field of the atomic nu- 
cleus and electrons) on N independent atoms. The sec- 
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ond one corresponds to the coherent radiation cross sec- 
tion. Let us denote (E ,P ), (Ex, Pi) and (i? 7 ,P 7 ) as 
the energy and momentum of incoming, outgoing elec- 
trons and emitted photon, respectively. The initial beam 
orientation with respect to the crystal axes is defined by 
two angles in the following manner. The three chosen or- 
thogonal axes of cubic crystal are [bi b2 hs\- The initial 
beam orientation is defined by the angle 9q between the 
initial electron momentum Po and crystal axis b3 and by 
the angle ipo between the electron momentum and crys- 
tal plane (b!b 3 ). Let 9 and tp be the emitted photon 
polar and azimuthal angles with respect to the direction 



of initial motion of the electron. Usually the polar angle 
is presented in the units of mc 2 /E in : 



E;, 



(2) 



where Ei n is the energy of initial particle, m is the elec- 
tron rest mass and c is speed of light. 

The angular-spectral distribution of CB (after integra- 
tion with respect to exit angles of electrons) is given by 
the following expressions [l|: 



J 



da 3 {x, O , «o, = ^rr^- —d£dtpl(x,6 ,a ,£„<p) 
2tt x 



(3) 



/(x, e , ao, c <p) = [i+(i- x) 2 ] (^i° h + *r c + - 2 (i - x) (v? h + *r + *f) 



(4) 



r 



where 7(x, 9q, ao, £, ip) is the intensity of radiation, N is 
the number of atoms in a crystal, er = Z 2 r^a, Z is the 
atomic number of medium, r is the classical electron 
radius, a is the fine structure constant, x = Ej/Eo is the 
relative energy of emitted photon and £ = 1/(1 + u 2 ). 

The i/'i'z.functions in equation Eq.|4]have the following 
structure III: 



*qjCoh 



coh 



g 

24^D s £(l-£)<7ifc 



(27r) 2 |5(g) | 2 [l-F(g)Y 



A 



Nn 



n 



eM-Ag 2 ) (5) 



where Sd is Dirac's delta function, g is the reciprocal 
lattice vector, A is the volume of an elementary cell of a 
direct lattice, 5(g) is the structure factor of the crystal, 
A is the root mean square of thermal displacement am- 
plitude of an atom from an equilibrium position. In our 
calculations we used Doyle- Turner parameterization Q 
for the atomic formfactor F(g). For a given orientation of 
the crystal, the longitudinal gn and transverse g± compo- 
nent of g with respect to the initial electron momentum 
Po are: 



where <?i,<?2, 33 are the projections of g on the crystal 
axes [bi bi 63] . The hS is a minimal value of the momen- 
tum transferred to medium along the direction of motion 
of primary particle: 



, hE^mc 2 
ho = 



(7) 



where h is the reduced Planck's constant. 



The incoherent 1PT2 functions in equation Eq. 2] have 
the following structure [![: 



y™ c = 6 + 4r(£) 

% nc = 6 + 24£(1 - £)r(0 



(8) 



where 



r(0 



1 1 mc 



/(C) +3(^/0 



(9) 
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33 cos 60 + (<7i cos ao + 32 sin ao) sin 9o 



(6) 



is the general expression for r(£) for arbitrary screening. 
The quantity $(H5/£) has the following form [|[: 



(q 2 - (M fie) 



dq 



(10) 
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where q is the momentum transferred to the nucleus. The 
quantity /(C) is the Coulomb correction: 



the atomic electrons, i.e. quantities tyf 2 were calculated 
by the theory given in Q . 



/(0 = c 2 E 



(ii) 



where C = Za. The corrections due to contributions by 



The threefold differential cross section of CPP by pho- 
tons in single crystals can be written in the following 
form [l[: 



Nag 
2tt 



dyd£d(pl(y,6a,aa,£,cp) 



(12) 



j( y , 6 Q , a , <p) = [y 2 + (i - y) 2 ] (n oh + *r + n l ) + 3 y(i - y) (*2 0/l + *r + 



(13) 



Here we use the following notation. The (i? 7 ,P 7 ) are 
the initial photon energy and momentum and (E + ,P + ), 
(£L,P_) are the energy and momentum of produced 
positron and electron pair. The y = E + / E-y is the rela- 



tive energy of produced positron, 



(E + /mc 2 )6 + and 



9 + is the angle between the positron and photon mo- 
menta. In case of CPP a minimal value of the momen- 
tum transferred to medium along the direction of motion 
of primary particle has the form: 



hE~,mcr 

no = 



2E+E- 



(14) 



The above described parameters are the same as for 
CB theory. 

Computed differential cross sections of CB (Eq. [3] in- 
tegrated over the emitted photon angles (6, ip)) and 
CPP (Eq. [12] integrated over the production angles (9± , 
ip±) of e + e~ pairs) are given in Fig. [TJ 

The top figure represents the spectral distribution of 
CB depending on the photon energy, x = Ej/Eq, on the 
Si and Ge single crystals. For an Eq = 178.2Gey electron 
beam making an angle of 9q — 30 mrad from the <001> 
crystallographic axis and ?po — 180/irad from the (110) 
plane of Si crystal, the maximum peak intensity occurs 
in the vicinity of 100 GeV as seen in Fig. [1] The angle 
ipo = 187 fir ad was chosen for the corresponding spectral 
distribution in case of Ge crystal. 

The differential cross sections of CPP in Si and Ge 
single crystals are given in the bottom of Fig. [T] depending 
on the energy of the e + e~ pair component, y = E±/E~ f . 
The total cross section of CPP has a maximum value 
for the chosen primary photon energy £L = 178. 2 GeV 
and orientation angles 6q = 30 mrad from the <001> 
crystallographic axis and ipo — O-dmrad from the (110) 
plane. 

Particles could undergo transmission at small angles 
with respect to the crystal axes and planes due to angular 
divergence and multiple scattering during the EMS de- 
velopment in crystal. CB and CPP processes have strong 
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FIG. 1: Differential cross sections of CB (top) as a function of 
photon energy, x = E^/Eq, and CPP (bottom) as a function 
of the energy of the e + e~ pair component, y = E±/E 1 for 
Si (solid curves) and Ge (dashed curves) single crystals. 



angular and energy dependences and the validity condi- 
tions of the Born approximation no longer hold at very 
high energies and small incidence angles with respect to 
the crystal axes and planes. The onset of this problem 
for the description of radiation emission and pair 
duction has the characteristic angle 8 V — Uo/mc 2 



4 



where Uq is the plane potential well depth. The CB and 
CPP theory may be applied for the incidence angles with 
respect to the crystal axes/planes 0q 3> V . The general 
theory of radiation and pair production is used [lol . |Tl| 
for incident angles 9q ~ @v and 9 < 9 V . 



III. MONTE CARLO SIMULATIONS 
A. General Considerations 

Monte Carlo simulation technique of EMS develop- 
ment in oriented single crystals is based on the coherent 
radiation and coherent e + e~ pair production processes at 
high energies. The program includes the formulae of co- 
herent and quasiclassical theories for radiation and pair 
production. The program calculated the cross sections 
directly using the formulae of coherent effects (Eqs. [3] 
and [T2"f . The quasiclassiacl theory was applied for calcu- 
lating the cross sections in case of small incident angles 
with respect to the crystal axes or planes. Direct calcula- 
tions by the formulae of this theory takes large computer 
time. Because the computer time needed to perform all 
the simulations is formidable, we used a data bank of pre- 
calculated cross sections of radiation and pair produc- 
tion to save the computer time. The data bank contains 
the total and differential cross sections on dependence of 
crystal type and its orientation, energy, horizontal and 
vertical angles and polarization. The cross sections for 
certain energy, angles, polarization were found by inter- 
polation of the numerical values of cross sections stored 
in a data bank. The interpolated values of cross sections 
differ from calculated values ~ 2 — 3%. We find that 
our approach yields very accurate numerical results for 
theory and a further saving of computer time. 



B. Simulation model 

The simulation model used in this work is based on 
random walk or particle history, where a history corre- 
sponds to following a particle from entering the medium, 
interacting and leaving it. Each particle history i can 
be represented by the array SJ [l2j denoting the state of 
particle before jth interaction: 



(r',Ar r r,).l-). q ') 



(15) 



where rl , O*- , rf- , E % - , q l - represent the electron or photon 
position, direction, polarization, energy and charge be- 
fore each interaction act, respectively. 

The initial particles, electrons or photons, generate two 
secondaries after each interaction. The initial electron 
loses energy via bremsstrahlung and produces secondary 
electron and bremsstrahlung photon. The electron can 
produce several photons during the traveling through 
matter. The number of produced photons depend on 
the thickness of the crystalline target and its orientation. 



The initial photon is transformed into the e + e~ pair af- 
ter interaction. The simulation code calculates new state 
of particles and state of produced particles after each in- 
teraction act. A history is terminated when the particle 
energy drops below a low energy cut-off, or when particle 
moves outside to the target. 

The similar simplified Monte Carlo simulation model 
was employed in [l3l fl4j . The simulation procedure 
has been improved by taking into account the peculiar- 
ities of the emission angles of photons and electrons in 
case of coherent effects (see Eq. [5]). This improvement 
gives the possibility to use the code for simulation of 
production of collimatcd photon beams from thin crys- 
tals (~ 20 — 100/im) [l5| . Polarization dependent cross 
sections and transverse dimensions of the particles beams 
are also implemented in the code. 



IV. SIMULATION RESULTS AND DISCUSSION 

A. Comparison with experiment 

A series of Monte Carlo simulations were performed 
for prediction of the results of CERN NA59 experiment. 
The goal of the experiment was the production of linearly 
polarized photon beams and conversion of the linear po- 
larization into circular with the help of oriented single 
crystals @, 1 M E3, EH • 
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FIG. 2: Total energy radiated for incoherent (top) and coher- 
ent (bottom) bremsstrahlung radiation in 1.5 cm thick silicon 
crystal pjLHq[. The solid curves represent Monte Carlo simula- 
tion results. The experimental results are: (■) - for unaligned 
silicon crystal, (A) - ^>o = 70 fir ad, (▼) - ipo — 180 /ir ad. 

The comparison of Monte Carlo simulation data with 
experimental data are given in Fig. O The top figure rep- 
resent the simulated and measured total radiated energy 
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TABLE I: Crystals and their orientations, radiation and ab- 
sorption lengths used in the simulations. 



Crystal EMS 



Orientation (mrad) Lr or La (cm) 



Si 
Si 
Ge 
Ge 



Initiated by e~ 
Initiated by 7 
Initiated by e~ 
Initiated by 7 



6>o=30; ^o=0.180 
00=30; Vo=0.6 
6*o=30; ^o=0.187 
0o=3O; Vo=0.6 



2.88 (9.36) 
5.45 (12.03) 
1.01 (2.30) 
1.76 (2.96) 



for 1.5 cm thick unaligned silicon crystal. The disoriented 
silicon crystal acts as an amorphous medium, hence we 
have incoherent bremsstrahlung spectra. 

The bottom figure represents the simulated and mea- 
sured CB spectra for the two orientations of silicon crys- 
tal with respect to incident electron beam. The elec- 
tron beam was incident at an angle of Qq = hmrad 
to axis <100> for both settings of silicon crystal. For 
the upper curve, the electron momentum makes an an- 
gle of tfjo = 70 firad with respect to the (110) plane and 
ipa = 180firad for lower curve. 

The experimental data are taken from 0, [Q. The 
prominent agreement between the results of simulation 
and experimental data is seen. The described Monte 
Carlo simulation code allows prediction of the photon 
spectra, linear and circular polarizations, optimal orien- 
tation and thickness of used crystals. 



B. Electromagnetic shower 

The developed code is capable to track the all gener- 
ations or histories of electrons and photons in crystal. 
Thus, it can be used for simulation of EMS develop- 
ment in single crystals oriented in CB mode. We con- 
sidered EMS initiated by the high energy electrons and 
photons in oriented silicon and germanium single crys- 
tals. The incident electron and photon beams parame- 
ters (energy, orientation angles, angular spread) are iden- 
tical with those in NA59 experiment and discussed above 
(Section [Hand Fig.E}. 

All calculations are carried out for the energy of ini- 
tial particles (electrons or photons) of 178.2 GeV. The 
simulations takes into account the initial electron beam 
angular divergence in both horizontal (48/j.rad) and ver- 
tical (33/irad) planes. The types of single crystals and 
their orientations used in Monte Carlo simulations are 
presented in Table HI A low energy cut-off of 5 GeV was 
placed on the production of all secondary particles. In 
this study no distinction is made between electrons and 
positrons (later on simply electrons). 

Important quantities in EMS development in matter 
are the radiation length (Lr) for charged particles and 
absorption length (La) for photons. There is a weak de- 
pendence of these quantities on particle energy in amor- 
phous media. Practically Lr and La are constant for a 
given amorphous material. There is a strong dependence 



of Lr and La on crystal type, its orientation, particle en- 
ergy and polarization in case of oriented single crystals. 
The corresponding radiation and absorption lengths for 
aligned and unaligned crystals (in parenthesis) are given 
in Table HI One can see a large reduction of L r and La 
in case of aligned crystals in comparison with amorphous 
media. 

Fig. [3] and Fig. Q] represent the Monte Carlo simula- 
tion results of longitudinal EMS development in silicon 
and germanium single crystals in energy range 5 GeV < 
E < 175 GeV . The energy carried by EMS is shown on 
the top figures in dependence on radiation and absorp- 
tion lengths. The bottom figures show the dependence of 
number of particles on radiation and absorption lengths. 
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FIG. 3: The energy (top) and number of particles (bottom) 
of EMS as a function of thickness of an aligned silicon sin- 
gle crystal. L r /a are the effective radiation or absorption 
lengths of the crystal for the EMS initiated by electrons or 
photons. Solid lines represents the EMS component initiated 
by electrons and dashed lines for EMS component initiated 
by photons. (•) - photon component of EMS, (■) - electron 
component of EMS. 



The energy carried by the photon component of EMS 
(initiated by electrons) reaches its maximum value at 
2 Lr for silicon and 1.6 Lr for germanium crystals. 
While the photon numbers reach the maximum at 3.5 Lr 
for silicon and 3.2 Lr for germanium crystals. The en- 
ergy carried by electrons is smoothly decreasing for both 
crystals. The number of electrons reach its maximum 
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FIG. 4: The energy (top) and number of particles (bottom) 
of EMS as a function of thickness of an aligned germanium 
single crystal. The notations are the same as for Fig. [3] 
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FIG. 5: Simulated photon (top) and electron (bottom) spec- 
tra, per unit of radiation length for a silicon crystal. The 
solid curves represent the spectra from aligned crystal and 
the dashed curves for unaligned crystal acted as amorphous. 



value at 3.5 Lr for silicon crystal and for germanium 
crystal at 4.2 Lr. The photon and electron components 
of EMS carry ~ 60% and ~ 30% of energy of initial elec- 
tron, respectively, for both crystals. The energy carried 
by the photon component of EMS are larger than the 
energy of electrons as seen from the figures. 

This is true for number of photons as well, which is 
about 2 times larger than the number of electrons. This 
is explained by the CB mechanism of radiation. Initial 
electrons lose large amount of energy around 100 GeV 
due to the CB radiation and produce high energy pho- 
tons as seen from Fig. [5] Simulated photon (top) and 
electron (bottom) spectra, Ed 2 N /cIEcILr, per unit of ra- 
diation length for a silicon crystal are shown. The solid 
lines represent the spectra from aligned crystal and the 
dashed lines for unaligned crystal acted as an amorphous. 
The photon spectrum for ICB smoothly decreases with 
increasing photon energy. One can observe an expressed 
photon peak in the vicinity of 100 GeV in case of CB 
spectrum. The initial electron loses energy mainly for 
producing high energy photons. Thus, the huge amount 
of energy is concentrated in the high energy region of 
the photon spectrum. This concentration of energy leads 
to the "delay" of EMS development. For example, the 
amorphous matter or crystal oriented in channeling set- 
tings split the energy of electrons very effectively and 



EMS develops more intensively ((J. 

The energy carried by the electron component of 
EMS (initiated by photons) reaches its maximum value 
at ~ 1.5 La for both crystals as seen from the Fig. [3]and 
Fig. |4] The number of electrons reaches the maximum 
value at ~ 3. La- The energy carried by the photon com- 
ponent of EMS are approximately the same as energy of 
electrons at their maximum values. But the number of 
photons is about 60% larger than number of electrons at 
La > 3. The behavior of the EMS initiated by photons 
is approximately the same as in the case of EMS initiated 
by electrons. The large number of electron-positron pairs 
are produced in the vicinity of 90 GeV due to the CPP. 



V. CONCLUSION 

The computer package presented in this paper is in- 
tended for Monte Carlo simulation of the electron and 
photon beams propagation through the oriented single 
crystals and amorphous media at high energies. The 
agreement between the simulation results and measure- 
ments is seen to be satisfactory. After finishing the pro- 
gram, the following parameters of the photons and elec- 
trons are reported: energy, angular and spatial distribu- 



7 



tions and polarization of passed electrons and photons. 
The important step in the code is the simulation of the 
photon radiation and e + e~ pair production angles, which 
strongly depends on the initial particles energy. This fea- 
ture of the computer code can be used for simulation of 
the photon beam collimation. 

The longitudinal behavior of electron and photon in- 
duced EMS in single crystals oriented in CB and CPP 
mode in the energy range between 5 GeV and 175 GeV 
was investigated. The interesting results are obtained 
concerning EMS development in single crystals aligned 
in CB or CPP mode. The crystal does not split the en- 
ergy of particles as amorphous media. It concentrates the 



energy of radiation and produced e + e~ pairs in the high 
energy region after the first interaction act. The crys- 
tals behave as a capacitor of energy up to a thickness of 
~ 2L R / A . As a result, the "delayed" EMS take place in 
the crystal. 
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